Phase transition in a one-dimensional Ising ferromagnet at zero-temperature under 
Glauber dynamics with a synchronous updating mode 
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In the past decade low-temperature Glauber dynamics for the one-dimensional Ising system has 
been several times observed experimentally and occurred to be one of the most important theoretical 
approaches in a field of molecular nanomagnets. On the other hand, it has been shown recently 
that Glauber dynamics with the Metropolis flipping probability for the zero-temperature Ising fer- 
romagnet under synchronous updating can lead surprisingly to the antiferromagnetic steady state. 
In this paper the generalized class of Glauber dynamics at zero-temperature will be considered and 
the relaxation into the ground state, after a quench from high temperature, will be investigated. 
Using Monte Carlo simulations and a mean field approach, discontinuous phase transition between 
ferromagnetic and antiferromagnetic phases for a one-dimensional ferromagnet will be shown. 

PACS numbers: 64.60.De Statistical mechanics of model systems, 64.60.-i General studies of phase transitions 
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I. INTRODUCTION 

Glauber dynamics for the Ising spin chain has been 
known for almost 50 years |lj , but only recently it became 
a really hot topic, not only from a fundamental, but also 
an applicative point of view [2 8]. It is well known that a 
purely one-dimensional (ID) system exhibits long-range 
ordering only at zero temperature T = OK . Nevertheless, 
in some situations long relaxation times for the magne- 
tization reversal with decreasing temperature can be ob- 
served, and finally at significantly low temperatures, the 
material can behave as a magnet. The phenomenon of 
slow magnetic relaxation is considered as one of the most 
important achievements of molecular magnetism, open- 
ing excitingnew perspectives including that of storing in- 
formation [9|, [Tfj| . Slow relaxation of the magnetization, 
predicted in the 1960s by Glauber in a chain of ferromag- 
netically coupled Ising spins [l|, in materials composed 
of magnetically isolated chains was observed for the first 
time in 2001 [2j. In 2002, this new class of nanomagnets 
was named single-chain magnets (SCM) [3j (for a recent 
review see Q) and the Glauber dynamics for the one- 
dimensional Ising spins system became one of the most 
important theoretical approaches for SCM. 

Within the Glauber dynamics for Ising spins with a 
spin s = 1/2, in a broad sense, each spin is flipped 
S%(t) — > — Si(t + 1) with a rate W(SE) per unit time 
and this rate is assumed to depend only on the energy 
difference implied in the flip. In this paper we consider 
the generalize class of zero-temperature dynamics defined 
as: 



W(SE) 
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which occurred to be very interesting not only from an 
applicative perspective, but also from a theoretical point 
of view as an example of no n-eq uilibrium dynamical sys- 
tems with many attractors [111 ] . The zero-temperature 
limits of the original Glauber dynamics [1| and Metropo- 
lis rates [l2[ (two the most popular choices) are respec- 
tively W G = 1/2 and W M = 1. 

Glauber dynamics was originally introduced as a se- 
quential updating (SU) process Also Monte Carlo 
method, used frequently for various models in statisti- 
cal physics, as proposed originally by Metropolis et al. 
[l^ . is essentially SU process. Evolution under dynam- 
ics defined by (fj} with random sequential updating is 
already well known in a case of one-dimensional sys- 
tem and can be derived analytically [111] . For any non- 
zero value of the rate Wq ferromagnetic steady state 
is reached and the dynamics belongs to the universal- 
ity class of the zero-temperature Glauber model [l|. 
The particular value Wo = corresponds to the con- 
strained zero-temperature Glauber dynamics ([ll| and 
references therein). In the constrained zero-temperature 
Glauber dynamics, the only possible moves are flips of 
isolated spins and therefore the system eventually reaches 
a blocked configuration, where there is no isolated spin 
fllj ] . i.e. for Wo = the relaxation time to the ferromag- 
netic steady state is infinite. 

The case of the synchronous updating, in which all 
units of the system are updated at the same time, is 
much more interesting. Moreover, clear evidence of a re- 
laxation mechanism which involves the simultaneous re- 
versal of spins have been shown experimentally for mag- 
netic chains at low temperatures [15|. 

In [2(| more general form of zero-temperature Glauber 
dynamics has been investigated than one defined by equa- 
tion by (TJ]). They have studied a model with two parame- 
ters r and S, which can be presented at T = analogously 
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where again W(5E) denotes the flipping rate per unit 
time. To fulfill the condition W(SE) S [0, 1], as seen from 
equation fSJ, the following relations have to be satisfied 



1 < 5 < 

r-2 



1 



r 

< J < 1 



(3) 



Above relations correspond to the region between thick 
lines in FigjTJ In 20] only the region denoted by the 
gray color in FigJJJhas been investigated (i.e. S < 0,T 6 
(0,1)). Comparing equations ^ and ([2]) we can easily 
derive the following relations: 
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In this paper we consider one-parameter model defined 
by Eq.(TT]) with W € [0,1]. Using relations @ we can 
determine corresponding regions in FigfT] (signatured by 
Wo € [0, 1]), which are disjoint from the gray region in- 
vestigated in [2(J. Also the area denoted by '?????' in 
FigfT] has not been investigated up till now - it could 
be considered using two-parameter model defined by ([2]) , 
but it is not covered by the one-parameter model, which 
is a subject of this paper. 



II. SIMULATION AND MEAN FIELD RESULTS 

We consider the chain of L Ising spins Oi = ±1 [i = 
1, 2, . . . L) with the periodic boundary conditions. In the 
initial state each lattice site is occupied independently by 
a randomly chosen value +1 or —1, both equally probable 
(high temperature situation) . In every time step all spins 
are considered simultaneously, but each spin is flipped 
independently with probability W(SE) defined by Eq. 
(UJ. It occurs that for all Wo S (0, 1) system eventually 
reaches one of the two final states - ferromagnetic steady 
state or antiferromagnetic limit cycle. If we measure the 
density of bonds (bond connects two sites with opposite 
spins) : 



1 L 



'i&i+l )■, 



(5) 



we obtain in the final state p s t — 1 (antiferromagnetic 
state) or p st — (ferromagnetic state). 

The time evolution of the mean value (averaged over 
10 4 samples) of the density of bonds measured in Monte 
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FIG. 1: Thick lines correspond to equations 5 — (1— T)/T and 
S = (r— 2) /r. The region between these two lines corresponds 
to the condition W(SE) € [0, 1]. In [2^| the region denoted by 
the gray color has been investigated (i.e. 8 < 0, F G (0, 1)), 
while in this paper we investigate two regions signatured by 
Wo 6 [0,1]. The region denoted by '?????' has not been 
investigated up till now, and it is not covered by the one- 
parameter model considered in this paper. 
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FIG. 2: The time evolution of the mean value of the density 
of bonds < p > measured in Monte Carlo steps for the lat- 
tice size L — 160 is presented. Averaging was done over 10 4 
samples. For Wo < 0.5 the mean number of bonds decreases 
in time to (ferromagnetic steady state) and for Wo > 0.5 
increases to 1 (antiferromagnetic limit cycle). 



Carlo steps (MCS) is presented in Figj2] This is seen 
that for Wo < 0.5 the average number of bonds decreases 
in time and eventually the system reaches the ferromag- 
netic steady state (< p(oo) >=< p st >= 0), while for 
Wo > 0.5 it increases and eventually antiferromagnetic 
limit cycle is reached (< p(oo) >=< p s t >= 1). Re- 
sults presented in Fig. [5] show that for Wo = 0.5 there is 
a phase transition between ferromagnetic and antiferro- 
magnetic phase. 

This phase transition can be predicted using the mean 
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field approximation analogously as it was done in [20j. 
In [20[ the mean field equations for the density of active 
bonds and magnetization has been derived: 



Bp + 2p 2 



= -5Tm(m z - 1). 

at 



(6) 



Using relations Q we can easily rewrite above equations 
in the case of our one-parameter model: 



^ = {l-2W )p{l-Bp + 2p 2 ) 

dm , 1 . , 
_ = (W - -)m(™ 2 
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As we see there are three types of fixed points: 

m st = ±1 and p st = 
m st = ±0 and p st = 1/2 
m st = ±0 and p st = 1 

In [20| only two first types have been considered: 

• p s t — (ferromagnetic state with m s t = — 1, 1) 

• Pst = 1/2 (so called active phase). 

However, there is a third fixed point p s t — l,w s t = 
0, which corresponds to antiferromagnetic steady state 
found in our computer simulations. It can be easily 
checked that: 

• for Wq < 0.5 ferromagnetic fixed point (m st — 
±l,p st = 0) is stable 

• for Wq > 0.5 antiferromagnetic fixed point (m st — 
0, p s t = 1) is the stable one. 

Thus, a mean field approximation predicts discontinu- 
ous phase transition between ferromagnetic and antifer- 
romagnetic phase for Wo = 0.5. This should be noticed 
that the transition value Wq = 0.5 corresponds to the 
original Glauber dynamics [1]. 

In the case of discontinuous phase transition one would 
expect the phase coexistence. We have provided com- 
puter simulations to confirm this mean field result and 
indeed coexistence of ferro- and antiferromagnetic phases 
can be observed near the transition point Wo = 0.5 (see 
FigEl). For Wo = 0.5 both types of clusters (ferro- and 
antiferromagnetic) are nearly the same size and after 
a long-time competition between them eventually one 
of two possible steady states is reached. Because for 
Wo = 0.5 both of them are equally probable we see the 
constant value of the average density of bonds in FigJ5] 
Let us now investigate the phase transition more quanti- 
tatively usin g M onte Carlo Simulations. 

Following |l4 , [20| we use as an order parameter the 
mean value of the density of bonds. We provide Monte 
Carlo simulations and wait until the system reaches the 
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FIG. 3: The time evolution of the Ising spins chain of the 
length L = 160 is presented. Black points represent bonds 
and thus black regions correspond to antiferromagnetic and 
white to ferromagnetic clusters. Coexistence of both types of 
clusters is visible for Wo w 0.5. For Wo = 0.5 both types 
of clusters are nearly the same size and there is a long-time 
competition between them leading eventually to one of two 
possible steady states (ferromagnetic or antiferromagnetic) 



final stationary state. Dependence between order param- 
eter in the stationary state < p s t > and the flipping prob- 
ability Wo is presented in Fig. [U showing again clearly 
discontinuous phase transition for Wq = 0.5 in agreement 
with the mean field result. In the case of Wo < 0.5 the 
ferromagnetic steady state is obtained with probability 
1 (for the infinite system L = oo ). For Wo > 0.5 the 
antiferromagnetic state is always reached, i.e. the sta- 
tionary states losses any remnants of the ferromagnetic 
Ising interactions. 

One of the most important issues connected with the 
coarsening is the relaxation time r, i.e. time needed to 
reach the ground state. In this paper we measure the re- 
laxation time starting from the random initial conditions 
and counting how many Monte Carlo steps is needed to 
reach the steady state (p = 1 or p = 0). We average 
over iV = 10 4 samples and calculate the mean relaxation 
time: 



1 



JV 



i=l 



(8) 



where Tj is the relaxation time of i-th sample. In Figj5] 
< r > divided by the square of the lattice size L as a 
function of the flipping probability Wo is shown. This is 
seen that for Wo = 0.5 the mean relaxation time scales 
as r ~ i 2 , which is well known result in a case of se- 
quential updating [l|| [ItJ • The dependence between the 
mean relaxation time < r > and the flipping probabil- 
ity Wo is non-monotonical. For Wo —> the relaxation 
time grows rapidly, which can be understood recalling 
that < r > if infinite for Wo = 11 1. For increas- 
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FIG. 4: Density of bonds pat in stationary state as a function 
of flipping probability Wo (so called exit probability) averaged 
over 10 4 samples. In the thermodynamical limit L — > oo for 
Wo < 0.5 ferromagnetic steady state is reached with proba- 
bility one (p st — 0) and Wo < 0.5 antiferromagnetic steady 
state is reached with probability one (p a t = 1)- Note that, the 
transition value Wo = 0.5 corresponds to the original Glauber 
dynamics. 
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FIG. 5: The mean relaxation times < r > divided by the 
square of lattice size L as a function of flipping probability 
Wo £ [0.48,0.52]. Averaging was done over 10 4 samples. Note 
that for Wo = 0.5 relaxation time scales with the system size 
as < t >~ L 2 . However, for Wo 7^ 0.5 scaling exponent 
differs from known value a = 2 (see Fig. 6). 



ing Wo the mean relaxation time decreases up to a cer- 
tain point W mm (L). However, due to the phase tran- 
sition in W = 0.5, for W 6 (W min (£),0.5) it grows 
again, resulting non-monotonic behavior shown in FigJS] 
The maximum peak is more and more narrow with the 
growing lattice size, which is expected behavior for the 
phase transition. The minimal value W mm (L) depends 
on the system size L as W™ m (L) = -2.5/ L + 0.5 and 
therefore lim£_ i . 00 W mm (L) — > 0.5. The mean relaxation 
time for this minimal value scales with the system size as 
< r(W mm ) >~ L 2 , i.e. with the same exponent as for 
the transition point Wo = 0.5. 

The most important question here is the one con- 
cerning the origin of the phase transition. As it was 
mentioned above, in the case of Metropolis flipping rate 
(Wo = 1) the system reaches antiferromagnetic limit cy- 
cle, instead for the ferromagnetic steady state [H|,[l4[. It 
can be easily understood, because for the flipping prob- 
ability Wo = 1, the case of synchronous updating is fully 
deterministic (see an example below): 

• • • mm • • • 

• • • mm • • • 

• • • mm • • • 

• • • mm • • • 

• • • mm ■■■ (9) 

On the other hand, only for Wo = 1 updating is really 
synchronous. For decreasing Wo only isolated spins are 
concerned really synchronously, since in the case of iso- 
lated spins SE < (see equation (fT])) the flip is provided 
with the probability 1. Flipping of isolated spins leads 



clearly to growth of ferromagnetic domains. Let us intro- 
duce for a while a notation L$e=o for the number of spins 
that flipping would not change the energy and Lse<o for 
the number of spins that flipping would decrease the en- 
ergy. The flip for SE = is realized with the probability 
Wo and for 5E < with the probability 1, which means 
that on average Lse<o + WqL$e=o is flipped in a single 
time step. In the case of Wo = 1, as mentioned above, the 
antiferromagnetic order is reached. On the other hand, 
for Wo = 1/Lse=o on average only one not isolated spin 
(i.e. with SE = 0) is flipped in a single time step, simi- 
larly to the case of the sequential updating for the system 
without isolated spins. Thus, because in the case of se- 
quential updating ferromagnetic steady state is reached, 
one can expect also ferromagnetic order in the case of 
synchronous updating for small values of W . Clearly 
the phase transition must occur somewhere between the 
antiferromagnetic order, preferred by a fully synchronous 
updating (Wo = 1), and the ferromagnetic steady state, 
preferred by sequential updating (Wo = \/L$e=o)- 

As mentioned above, for Wo = 0.5 and Wo = W rrim (i) 
the mean relaxation time scales with a system size as 
~ L 2 . We have checked also the scaling for other values 
of Wo and we have obtained power laws < r(Wo) >^ L a 
with Wo-depending scaling exponents a — a(Wo). The 
dependence between scaling exponent and the flipping 
probability is presented in Fig. [5] The shape of the 
curve a(Wo) mimic the shape of < r(Wo) >, which can 
be understood looking at Fig. O 
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FIG. 6: The mean relaxation time < r > scales with the 
system size as < r >~ L a . For Wo = 0.5 the scaling exponent 
a — 2, which is well known result in the case of sequential 
updating. However, in general scaling exponent depends on 
the flipping probability Wo, i.e. a — a(Wo). Dependence 
between scaling exponent a and the flipping probability Wo 
is shown. Simulations were done for the system size L £ 
[20, 1280] and averaged over 10 4 samples 



III. SUMMARY 

In this paper we have been investigating the relaxation 
of the Ising spins chain under the generalized class of 
Glauber dynamics at zero-temperature. Within such a 
dynamics, the flipping probability in a case of conserved 
energy is given by arbitrary value of Wo € [ 0, 1 ] (review in 
a case of sequential updating can be find in ) . We have 
proposed to use synchronous updating for such a gener- 
alized class of zero-temperature dynamics. Our motiva- 
tion for this work came from recent experiments showing 
slow relaxation in magnetic chains at low temperatures 
0-H, [HI- We have shown by Monte Carlo simulations 
that there is a phase transition for Wo = 0.5, which cor- 
respond to the value originally proposed by Glauber 
Following [2(| we were able to obtain the mean field result 
which predicts discontinuous transition between ferro- 
and antiferromagnetic phases for Wq = 0.5. 
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